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Abstract A massive star is defined to be one with a mass greater than ~ 8 — 10A4 Q . 
Central to the on-going debate on how these objects [massive stars] come into being is 
the so-called Radiation Problem. For nearly forty years now, it has been argued that the 
radiation field emanating from massive stars is high enough to cause a global reversal of 
direct radial in-fall of material onto the nascent star. We argue that only in the case of 
a non-spinning isolated star does the gravitational field of the nascent star overcome the 
radiation field. An isolated non-spinning star is a non-spinning star without any circum- 
stellar material around it, the gravitational field beyond its surface is described exactly 



CZ2 ' by Newton's inverse square law. The fact that massive stars should have their gravita- 



tional field being much stronger than their radiation field is drawn from the analysis of 



an isolated massive star, where in this case the gravitational field is much stronger than 



the radiation field. This conclusion is erroneously extended to the case of massive stars 
enshrouded in gas & dust. We find that, for the case of a non-spinning gravitating body 
where we take into consideration the circumstellar material, that at ~ 8 — 10.M©, the 
radiation field will not reverse the radial in-fall of matter but a stalemate between the radi- 
ation and gravitational field will be achieved - i.e., in-fall is halted but not reversed. This 
picture is very different from the common picture that is projected and accepted in the 
Z^ popular literature that at ~ 8 — 10A4 ©, all the circumstellar material - from the surface 

f~*l , of the star, right up to the edge of the core; is expected to be swept away by the radiation 

£f) • field. We argue that massive stars should be able to come into being if the molecular core 

from which they form exhibit some rotation because a rotating core exhibits an ASGF 
f—j ■ which causes there to be an equatorial accretion disk and along this equatorial disk, the 

radiation field can not be much stronger than the gravitational field hence this equatorial 
accretion disk becomes the channel via which the nascent massive star accretes all of its 
material. 



Key words: (stars:) circumstellar matter - (stars:) formation - radiative transfer. 
PACS (2010): 97.10.Bt, 97.10.Gz, 97.10.Fy 



1 INTRODUCTION 

According to current and prevailing wisdom, it is bona-fide scientific knowledge that our current under- 
standing of massive star formation is lacking and this is due to the existing theoretical and observational 
dichotomy. In the gestation period of a star's life, its mass will grow via the in-falling envelope (i.e., 
circumstellar material) and also through the forming accretion disk being along it's equator. As far as 
our theoretical understanding is concerned, this works well for stars less than about 8 — 10.M©. In the 
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literature, it is said that the problem of massive stars (M. s tar > 8 — lCM©) arises because as the central 
prostar's mass grows, so does the radiation pressure from it, and at about 8—10 A4q, the star's radia- 
tion pressure becomes powerful enough to halt any further in-fall of matter onto the protostar (Larson 
& Sta rfield [1 9711 Kahn[197l Yorke<& KrugelQjjfl Wolfire & Cassinelli [19871 Palla & Stahler[1993 
Yorke [200121 Yorke & Sonnhalter |2003i l. So the problem is; how does the star continue to accumulate 
more mass beyond the 8 — 10 A4q limit? If the radiation field really did reverse any further in-fall of 
matter and protostars exclusively accumulated mass via direct radial in-fall of matter onto the nascent 
star and also via the accretion disk, this would set a mass upper limit of 8 — 10 Mq for any star in 
the Universe. Unfortunately (or maybe fortunately) this is not what we observe. It therefore means that 
some process(es) responsible for the formation of stars beyond the 8—10 M.q limit must be at work, 
hence, a solution to the problem must be sought because observations dictate that it exists. 

If this is the case, i.e., the radiation problem really did exist as stated above, and our physics where 
complete viz gravitation and radiation transport, then, the solution to the conundrum would be to seek 
a star formation model that overcomes the radiation pressure problem and at the sametime allowing for 
the star to form (accumulate all of its mass) before it exhausts its nuclear fuel. Two such (competing) 
models have been set-forth, i.e., (1) the Accelerated Accretion Model (AAM) (Yorke [20021 12003b and, 
(2) the Coalescence Model (CM) (Bonnell et al. 11^511211021120061120071 

The latter scenario, i.e. the CM; is born out of the observational fact that massive stars are generally 
found in the centres of dense clusters (see e.g. Hillenbrand 119971 Clarke et al. 120001 ). In these dense 
environments, the probability of collision of proto-stellar objects is significant, hence the CM. This 
model easily by-passes the radiation pressure problem and despite the fact that not a single observation 
to date has confirmed it (directly or indirectly), it [CM] appear^] to be the most natural mechanism by 
which massive stars form given the said observational fact about massive stars and their preferential 
environment. 

The AAM is just a scaled up version of the accepted accretion paradigm applicable to Low Mass 
Stars (LMSs). This accretion takes place via the accretion disk and for the reason mentioned above 
that the accretion mechanism must be such that it allows for the star to form before it exhausts its 
nuclear fuel, the accretion can not take place at the same steady rate as in the case of LMSs [A4 < 
3Mq) but must be accelerated and significantly much higher. While there exists many examples of 
massive stars surrounded by accretion disks, one of the chief obstacles in verifying this paradigm is that 
examples of HMSs tend to be relatively distant (> 1 kpc), deeply embedded, and confused with other 
emission sources (see e.g. Mathews et al. 120071) . Additionally, HMSs evolve rapidly, and by the time an 
unobstructed view of the young star emerges, the disk and outflow structures may have been destroyed, 
consequently, observations to date have been unable to probe the 10 — 100 AU spatial scales over which 
outflows from the accretion disks are expected to be launched and collimated (e.g. Mathews et al. 12007b . 

The other alternative, which is less pursued, would be to seek a physical mechanism that overcomes 
the radiation pressure problem as has been conducted by the authors Krumholz et al. ( 2005 2009). These 
authors (Krumholz et al. 2005 2009 ) believe that the radiation problem does not exist because radiation- 
driven bubbles that block accreting gas are subject to Rayleigh-Taylor instability which occurs anytime 
a dense, heavy fluid is being accelerated by light fluid for example when a cloud receives a shock, or 
when a fluid of a certain density floats above a fluid of lesser density, such as dense oil floating on water. 
The Rayleigh-Taylor instabilities allows fingers of dense gas to break into the evacuated bubbles and 
reach the stellar surface while in addition, outflows from massive stars create optically thin cavities in 
the accreting envelope. These channel radiation away from the bulk of the gas and reduce the radiation 
pressure it experiences. In this case, the radiation pressure feedback is not the dominant factor in setting 
the final size of massive stars and accretion will proceed albeit at much higher rates. Amongst others, 
the model by the authors Krumholz et al. ([2005, 2009) is mechanical rather than natural, in that Nature 
has to make a special arrangement or must configure herself in such a way that massive stars have a way 
to come into being. Does there not exist a smooth and natural way to bring forth massive stars into this 
World? 

1 This is on the assumption that our understanding of gravitation and radiation transport is complete. 
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In this reading, we redefine the radiation problem (for the spherically symmetric case) and we do 
this via a subtle and overlooked assumption made in the analysis leading to the radiation problem; that 
the surroundings of the protostar is a vacuum (see e.g. Yorke 2002. Yorke & Sonnhalter 2003; Zinnecker 
& Yorke [2002b ; surely, this is clearly not true. The researchers Yorke [20021 Yorke & Sonnhalter [20021 
Zinnecker & Yorke[20O7; amongst others, hold the view that; from a theoretical stand-point, the radiation 
field is stronger than the gravitational field for massive stars hence the in-fall process of material must 
be reversed; but this conclusion has been reached - as will be shown in the next section; after comparing 
the gravitational field strength at point r of a star in empty space to its radiation field strength at point r. 
In practice, stars are found embedded inside a significant mass of gas and dust. The radiation problem 
is arguably the most important problem of all in the study of the formation of stars, hence thus, it is 
important to make sure that this problem is clearly defined and understood. 

Having taken into consideration the circumstellar material, we find that at ~ 8 — lOA^©, the radi- 
ation field will not reverse the radial in-fall of matter but a stalemate between the radiation and gravi- 
tational field will be achieved - i.e., in-fall is halted but not reversed. Certainly, this picture is not at all 
congruent (or somewhere near there) to the common picture that is projected and accepted in the popular 
literature where at ~ 8 — IO.M0, all the circumstellar material - from the surface of the star, right up 
to the edge of the core; is expected to be swept away by the powerful radiation field. This finding is 
not a complete but a partial solution to the radiation problem in that beyond the 8 — 10.M© limit, the 
nascent star will not accrete any further, its mass will stay put at this value each time accreting from the 
stagnant and frozen envelope once its mass drops below this 8 — IOA^q limit. Very important to note 
is that this is for a spherically symmetric gravitational setting where the gravitational field has only the 
radial dependence and is described exactly by Newton's inverse square law. 

In a different reading i.e., Nyambuya ( l2010aD an Azimuthally Symmetric Theory of Gravitation 
(ASTG) was set-up and there-in a thesis was set-forth to the effect that; (1) for a non-spinning star, 
its gravitational field is spherically symmetric, i.e., it is only dependent on the radial distance from 
the central body; (2) for a spinning gravitating body, the gravitational field of the body in question is 
azimuthally symmetric, i.e., it is dependent on the radial distance from the central body and as-well 
the azimuthal angle. In a follow-up reading i.e., Nyambuya (1201061 l: it has been shown that the ASTG 
predicts (1) that bipolar outflows may very well be a purely gravitational phenomenon and also that; 
(2) along the spin-equator of a spinning gravitating body, gravity will channel matter onto the spinning 
nascent star via the spin-equatorial disk without radiation having to reverse this inflow, thus allowing 
stars beyond the critical mass 8 — 10A4 0, to come into being. 

If the ASTG proves itself, then the present reading together with Nyambuya (2010a 20106) com- 
prise (in our view) a solution to the radiation problem. Given that the solution to this problem has been 
sought via sophisticated computer simulations and lengthy numerical solutions, and additionally, given 
the simplicity and nai've-ness of the present approach which seeks to further our understanding of this 
problem; perhaps - this reading presents not only my misunderstanding of the problem, but also of the 
approach to the problem - but more on the optimistic side of things, I believe the radiation problem has 
here been understood and that the approach is mathematically and physically legitimate, so much that 
we are of the objective view that to they that seek a solution to this problem, this reading is something 
worthwhile. 



2 THE RADIATION PROBLEM 

Following Yorke [20021 for direct radial accretion and accretion via the disk to occur onto the nascent 
star, explicitly, it is required that the Newtonian gravitational force, GA4 s tar(t)/r 2 , at a point distance 
r from the star of mass Ai s tar and luminosity C star {t) at any time t, must exceeds the radiation force 

K e ffC sta r{t)/^cr 2 i.e.: 

GA4 s tar(t) Keff£star(t) 

r 2 Ancr 2 
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where c = 2.99792458 x 10 8 ms _1 is the speed of light in vacuum, G = 6.667 x 1CT 11 kg _1 m 3 s -2 
is Newton's universal constant of gravitation, K e ff is the effective opacity which is the measure of the 
gas's state of being opaque or a measure of the gas imperviousness to the rays of light and is measured 
in m 2 kg~ 1 . This analysis by Yorke d20021 > which is also reproduced in Zinnecker & Yorke d2007l l, is a 
standard and well accepted analysis that assumes spherical symmetry and at the sametime it does not 
take into account the nascent star's circumstellar material. On the other hand, star formation is not a 
truly spherically phenomenon (see e.g. reviews by Zinnecker & Yorke |2"0"0"T1 McKee & Ostrikker l2007ft 
but this simple calculation suffices in as far as defining curtain-region of 8 — lOA^© when radiation 
pressure is expected to become a significant player on the star formation podium. What will be done in 
this reading is simple to perform the same calculation albeit with the circumstellar material taken into 
account. In the penultimate of this section, we shall make our case based on the just said. 

Now, this calculation by Yorke (12002b and Zinnecker & Yorke (12007b . proceeds as fol- 
lows; the inequality dTJ, sets a maximum condition for accretion of material, namely K e ff < 
4ir cG Ai s tar (t)/£ s tar (t), and evaluating this we get: 

Ke// <1.30xl0^(^gW)(^)" 1 , (2) 

where M. s tar{t) and C s tar{t) are in solar units. Given that, £ star (t) = £ Q (M. s ta r (t) / Mq) , this 
implies that: 

/1qn 1 n4 ( M star (t) \ ^ ( Mstar \ ^ ( 00X104 \ 1/2 

K e ft < 1.30 x 10 — , , =$■ —rz — > ■ (3) 

" \ Mo J \ M Q J \ K e ff J 

Now, given that the dusty Interstellar Medium's (ISM) averaged opacity is measured to be about 
20.0 m 2 kg _1 (Yorke l2~002b and using this (as an estimate to setting the minimum critical mass, see 
Yorke [20021 Zinnecker & Yorke [2"007l l. we find that this sets a minimum upper mass limit for stars of 
about 10.M© for gravitation to dominate the scene before radiation does. It is clear here that the opacity 
of the molecular cloud material is what sets the critical mass, thus a cloud of lower opacity will have 
a higher critical mass. It is expected that the opacity inside the cloud will be lower than in ISM, hence 
thus, in adopting the value K e ff = 20.0m 2 kg _1 (see Yorke [20021 Zinnecker & Yorke [2007l l. this was 
done only to set a minimum lower bound for massive stars. Dust and gas opacities are significantly 
frequency-dependent and one has to take this into account for a more rigid setting up of a minimum 
mass for when the radiation field is expected to overcome the gravitational field. 

As can be found in Yorke (12002b . the AAM finds some of its ground around the alteration of the 
opacity. For example, if the opacity inside the gas cloud is significantly lower than the ISM value, then 
accretion can proceed via the AAM. To reduce the opacity inside the gas & dust cloud, the AAM posits 
as one of the its options that optical and Ultra- Violet (UV) radiation inside the accreting material is 
shifted from the optical/UV into the far Infrared (IR) and also that the opacity may be lower than the 
ISM value because the opacity will be reduced by the accretion of optically thick material in the blobs of 
the accretion disk. Thus reducing the opacity or finding a physical mechanism that reduces the opacity 
to values lower than the ISM is a viable solution to the radiation problem. The above mechanism to 
reduce the opacity are rather mechanical and dependent on the environment. 

Now that we have presented the radiation problem as it is commonly understood, we are ready to 
make our case by inspecting (Q]). Clearly and without any doubt, the left hand side of this inequality 
is the gravitational field intensity for a gravitating body in empty space while the right hand side is 
the radiation field of this same star in empty space. From this - clearly; we are actually comparing the 
radiation and gravitational field intensity of a star in empty space, whereas the real setting in Nature, 
stars are found heavily enshrouded by gas and dust. Clearly, the conclusions that one finds from (Q3 such 
as that - at about 8 — ICMq the radiation field of the nascent star is powerful enough to not only halt but 
reverse the in-fall of material onto the nascent star; this can not be extended to the scenario where a star is 
submerged in gas and dust, it is erroneous to do so. Clearly, at this very simplistic, naive and fundamental 
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level, there is a need to redefine the radiation problem by including in the left hand side of (Q]l, the 
circumstellar material. Wolfire & Cassinelli (119871) amongst others, have performed this calculation 
where they have taken into account the circumstellar material and reaching similar conclusions (as e.g. 
those of Yorke 2002]i. We reach a different conclusion to that of Wolfire & Cassinelli (1987) because 
unlike these researchers; we use the observational fact that molecular clouds and molecular cores are 
found exhibiting a well behaved density profile p oc r~ a <>, and from this, we calculate from this a 
general mass distribution (Ai oc r~ a ) and we use this to compare the gravitational and radiation field 
strengths at point r and from there draw our interesting conclusions. 



3 RADIATION AND THE CIRCUMSTELLAR MATERIAL 

Neglecting thermal, magnetic, turbulence and any other forces (as will be shown latter on in this section, 
these forces do not change the essence of our argument, hence there is no need to worry about them here) 
and considering only the gravitational and radiation field from the nascent star, we assume here that a 
star is formed from a gravitationally bound system of material enclosed in a volume space of radius 
TZ-core{t) and we shall call this system of material the core and further assume that this core shall have a 
total constant mass M. cor e at all times. Now for as long as the material enclosed in the sphere of radius 
r < lZ C ore(t) is such that: 



GM(r,t) K eff £ st ar(t) 



> 



47TCX 2 



(4) 



then, radiation pressure will not exceed the gravitational force in the region r < lZ core (t) hence thus 
direct radial in-fall is expected to continue in that region. If M cs i(r, t) is the mass of the circumstellar 
material at time t enclosed in the region stretching from the surface of the star to the radius r, then, 
M.(r,t) = A4 cs i(r,t) + M s tar(t), hence the difference between Q and (Q]) is that in (|4|i we have 
included the circumstellar material. This is not the whole story. 
Now, (O can be written differently as: 



M(r,t) > 



K 



ef f t^star 



(f) 



4ttGc 



(5) 



which basically says as long as the amount of matter enclosed in the region of sphere radius r satisfies 
the above condition, the radiation force will not exceed the gravitational force in that region of radius r. 
In fact, (|5]l is the Eddington limit applied to the region of radius r. This is identical to equation (10) in 
Wolfire & Cassinelli (119871 1. In their work, Wolfire & Cassinelli dl987l l solve numerically the radiative 
transfer problem to determine the effective opacity at the outer edge of the massive star forming core and 
from this they determine the limits for the grain-sizes that are needed for the formation of massive stars. 
Wolfire & Cassinelli d!987b 's approach is a typical approach used to probe the conditions necessary for 
massive stars to form. 

Our approach is very different from that of Wolfire & Cassinelli (4987) and most typical approaches 
used to study the radiation problem where sophisticated computer simulations and numerical solutions 
are used. Ours is a simple and naive approach needing no computer simulations nor numerical codes. 
We shall insert A4(r, t) = Ai cs i(r, t) + M. 8 tar(t) into (|4|i and thereafter rearranging, one obtains: 



Mcsi{r,t) > 



H 



ef f *^star 



(t) 



4ircGM st ar(t) 



- 1 



M st ar{t) = 



Mstarjt) 
IOMq 



- 1 



M st ar(t), 



(6) 



and our main thrust is to seek values of r in the above inequality that satisfy it. We shall do this by 
finding a form for A4 cs i (r, t). 

Before doing this, let us apply © to the entire core i.e., r — lZ CO re- This must give us the condition 
when the star's radiation field is strong enough to sweep away all the circumstellar material from the 
surface of the star right up to the end of the core; so doing, one finds that the star's luminosity should be 
such that: 
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K e ff£ atar (t) 

AnGc 

In making this calculation, we have made the tacit and fundamental assumption that the star's mass will 
continue to increase until the star reaches a critical luminosity determined by the mass of the core - let 
us denote this critical luminosity by C* core . From the above, it follows that: 

4ncGM CO re „. 

'-'core ~ (°) 

K eff 

With this defined, then for the radiation field to globally overcome the gravitational field, the nascent 
star's luminosity must exceed the critical luminosity of the core, i.e.; 

C star (t) > £* core . (9) 

Now, knowing the mass-luminosity relationship of stars is given by £ sta r(t) = £© (M (t)/M © ) , 



then the critical condition C s t a r (t) = C* core will occur when 

'M st ar\ _ ( Keff£ & \ "^ ( M 



V3 / kA X 1/3 

^1 . (10) 

Given this and taking K e ff = 20.0 m 2 kg _1 and then plucking this and the other relevant values - G, c, 
etc; into the above, we are lead to: 

M c ^ /3 



/ JW-rrbax \ \ " ' * i:ui ( - 1 



1 core. 



V M Cl) ) ViO-^c 

where we have set A4 s tar = -M-max- As already said, using K e ff = 20.0 m 2 kg _1 gives us the minimum 
lower bound. What this means is that the mass of the core from which a star is formed may very well be 
crucial if not critical or detrimental in deciding the final mass of the star because the mass of the core 
determines the time when global in-fall reversal will occur. 

From this simplistic and rather naive calculation, we can estimate the efficiency of the core: 

thus a 100.Mo core will (according to the above) form a star at an efficiency rate of about 2% and it will 
produce a star of mass 2_A/f©. A 10A1© star will be produced by a core of mass lO 4 ^© at an efficiency 
rate of about 0.1%. A 10 4 .M © core is basically a fully-fledged molecular cloud. The production of this 
IO.M0 star is on the assumption that the rest of the material (i.e., 10 4 .Mq - 10A^© = 9.99 x 10 4 Mq) 
will not form stars. In reality, some of the material in this 10 4 A^© core will form many other stars. 
Further, a lOOA'f© star will form in a GMC of mass about 10 7 .M©. The above deductions that high 
mass stars will need to form in clouds of mass > 10 4 7W©, resonates with the observational fact that 
massive stars are not found in isolation (e.g. Hillenbrand 119971 Clarke et al. 120001) since the other 
material will form stars. 

The relationship ( fTTT i is interesting viz its similarity to Larson's 1982 empirical discovery. With a 
handful of data, Larson (1982 ) was the first to note that the maximum stellar mass of a given population 
of stars is related to the total mass of the parent cloud from which the stellar population has been born. 
That is to say, if M c i is the mass of molecular cloud and M max is the maximum stellar mass of the 
population, then: 

m -=(^T (i3) 
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where M.q — 13.2A4© and ctL = 0.430. This law was obtained from a sample of molecular clouds 
whose masses are in the range 1.30 < log 10 (M./ M.q) < 5.50. Larson's Law is thought of as being a 
result of statistical sampling but we are not persuaded to think that this is the case, such a coincidence 
is - in our opinion and understanding, to good to be true. We believe Larson's Law is Nature's subtle 
message to researchers; it is telling us something about the underlaying dynamics of star formation. 
This said, could the relationship ( fTTT ) be related to Larson's result? The indices of Larson's relation and 
relationship ( fTTT i have a deviation of about 33% and the constant Mo has a similar deviation of about 
33%. Could Larson's fitting procedure be "tuned" to conform to relationship ( fTTT ) and if so, does that 
mean Larson's relationship finds an explanation from this? 

Perhaps the deviation of our relation from that of Larson may well be that our result is derived from 
an ideal situation where we have considered not the other forces such as the magnetic, thermal forces 
etc, also, we have considered star formation as a spherically symmetric process of which it is not and this 
may also be a source of correction to this result in order to bring it to Larson's result. Let us represent 
all these other forces by Mother (<?.g. magnetic, turbulence, viscous etc). Clearly these forces will not 
aid gravity in its endeavor to squeeze all the material to a single point but aid the radiation pressure in 
opposing this. Given this, it means we must write (0]i as: 

GM(r,t) K e ffC st ar(t) [Fptfeerj .... 

r 2 Ancr 2 m 

where m is the average mass of the molecular species of the material constituting the cloud. The above 
can be written in the form: 

- M ^ncG{M(r,t)-r*\F other \/m) 

K eff 

and writing Ai'(r, t) = r 2 |F ot / ler |/rn, we have from the above: 

^ncG[M(r,t)-M'(r,t)] 

L s tar{t) < , (16) 

K eff 

and from this it is clear that the other forces will act in a manner as to reduce the critical luminosity of 
the core thus our result (fTTb . when compared to natural reality where these other forces are present, it is 
expected that a deviation from the real observations must occur. As stated in the opening of this section 
that the inclusion of the magnetic, thermal forces etc will not change the essence of our argument, hence 
the above justifies why we did not have to worry about these other forces as the essence of our result 
stands. The situation is only critical when these other forces become significant in comparison to the 
gravitational force. 

In the succeeding section, we compute the mass distribution function and from there show that 
one arrives at the same result as ©. Additionally and more importantly, we are able to compute the 
boundaries where the radiation field will be strong enough to overcome the gravitational field. Amongst 
other interesting outcomes, we shall see that the radiation field will create a cavity inside the star forming 
core and that this cavity grows with time in proportion to the radiation field of the nascent star. 

4 MASS DISTRIBUTION FUNCTION 

First we compute the enclosed mass M.{r,t). We know that stellar systems such as molecular clouds 
and core are found exhibiting a radial density profiles given by: 

p{r,t)=p (t)( r -^X P (17) 



where po (i) and tq (t) are time dependent normalization constants and a p is the density index. In order 
to make sense of this density profile ( fT7T > we shall have to calculate these normalization constants. In 
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its bare form, the power law ( fT71 ) as it stands implies an infinite density at r = 0. In general, power 
laws have this property. Obviously, one has to deal with this. The usual or typical way is to impose a 
minimum value for r, say r = r m ,;„ = ro(i) and, assign a density there. Here, this minimum radius has 
been made time dependent for the sole reason that if the cloud is undergoing free fall as in the case in 
star formation regions, this quantity will respond dynamically to this, hence it will be time dependent. 
Now, for a radially dependent density profile, the mass distribution is calculated from the integral: 



M(r,t)= 4irr 2 p(r,t)dr. (18) 

Inserting the density function (TTTI i into the above integral and then evaluating the resultant integral, we 
are lead to: 

M(r,t) = ( ^f^ ) (r*-' - r^(*)) , (19) 

and this formula does not apply to the case a p = 3 hence this is valid for < a p < 3. The case 
a p = 3 is described by a special MDF which is: A4(r, t) = [47rpo(£) r mmW] m ( r / r mm(t))- We shall 
not consider this case as it will not change the essence of our argument. 

Now, what we shall do here is to constrain a p and show that: < a p < 3. This exercise is being 
conducted to define the domain which our result has physical significance. First we shall establish that 
Up < 3 and this we shall do by using the method of proof by contradiction. Let (n > ri). For this 
setting, we expect that [M.(r 2 ) > A4(ri)] and this is obvious thing because as we zoom out of the 
molecular cloud radially, one would expect to have in the bigger sphere of radius r 2 more matter than 
that enclosed in the smaller sphere of radius n; therefore our conditions is: [A4(r 2 ) > A4(ri) => 
M(r 2 ) — M(ri) > 0]. Using equation ( fT9l l. we have: 

M(r 2 ) - M(n) = ( An /° r ™A (V 2 3 ~^ _ r 3-M > 0, (20) 



andfor(a p > 3) we have (3— a p < 0)so when we divide by the term (47rpQr^ lin )/(3—a p ) on bothsides 
of the inequality, we must change the sign of the inequality from > to < because (47rpo?' m 'i„)/(3 — a p ) 
is a negative number. So doing, we will have: r 2 ap —r x ap < 0, and this implies r 2 Qp < r x ap and 
from this follows directly the relationship: 

ri < r 2 , (21) 

and this is a clear contradiction because it violates our initial condition [r 2 > r\ => M(r 2 ) > M. (ri)] 
as this is saying [r 2 < r\ ==> Ai(r 2 ) > M(ri)] which is certainly wrong. From a purely mathematical 
stand-point, we are therefore forced to conclude that a p < 3 if the condition [r 2 > r\ => Ai(r 2 ) > 
M(n)] is to hold -QED. 

Now we shall establish that a p > and we shall do this using physical arguments. If 3 — a p > 3, it 
means as one zooms out of the cloud from the center, the cloud's average material density increases. This 
scenario is unphysical because gravity is an attractive inverse distance law and thus will always pack 
more and more material in the center than in the outer regions hence thus the only material configuration 
that can emerge from this setting is one in which the average density of material decreases as one zooms 
out of the cloud. This implies 3 — a p < 3 which leads to a p > 0, hence combining the two results, we 
are going to have: < a p < 3. Now we have defined the physical boundaries of the density profile. 

Now we have to normalize the MDF by imposing some boundary conditions. The usual or tradi- 
tional boundary condition is to set Ai(r m i n , i) = and this in actual fact means there will be a cavity 
of radius r m i n (t) in the cloud. What we shall do is different from this normal or traditional normal- 
ization. We shall set M.(r m i n ,t) = M. s tar where M. s tar is the mass of the central star hence thus 
fmin(t) = T^-star{t). Thus what we have done is to place the nascent star in the cavity. This means we 
must write our MDF as: 
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= ^P^XtaM } / 3-a„ -TllZ"{t 



and this applies for TZ s tar{t) < r < T^core(t). 

Now, if the mass enclosed inside the core remains constant throughout, then we must have at 
r = 7lcore(t) the boundary condition M.(TZ C ore, t) = Mcore- We know that the sum total of all the 
circumstellar material at any given time is given by: M C siif) = -A4 core — M-star{t)- Combining all the 
information, we will have: 

n7T P0 (t)C(t) \ = Mcsljt) 

and this means the MDF can now be written as: 



Circumstellar Material in Region Radius r 



Mass of the nascent star 



M{r,t) = Mcsi(t) 3 _ a /n sta 3 r _ a v - + M s tar(t) for r > Tl star (t). (24) 

We shall take this as the final form of our mass distribution function. If the reader accepts this, then 
what follows is straight forward exercise and leads to what we believe is a significant step forward in 
the resolution of the radiation problem. The reader may want to query that we have overstretched our 
boundary limits by making the MDF be continuous from the surface of the star right up till the edge of 
the core. In that event, we need to clear this and reach an accord. 

First let us consider a serene molecular core way before a star begins to form at the centre. We 
know that the density is not a fundamental physical quantity but a physical quantity derived from two 
fundamental physical quantities which are mass and volume, i.e., density=mass/volume: we must note 
that this is defined for (volume> 0). We shall assume that this core exhibits the density profile p oc r~ a f> . 
This fact that p oc r~ ap combined with the fact that density not a fundamental physical quantity but 
quantity derived from two fundamental quantities, suggests thats at any given time the mass must be 
distributed in proportion to the radius, i.e., A4(r, t) oc r a . The radial dependency of the density is an 
indicator that that mass has a radial dependency. The relationship A4 (r, t) oc r a means we must have 
M.(r, t) — ar a + b where (a, b) are constants. We expect that M.(Q, t) = 0. If this is to hold (as it must), 
then (b = 0) and (a > 0). We also expect the condition M.(1Z core , t) — M. CO re to hold. If this is to hold 
(as it must), then we will have a = M core /TZ" ore (t) hence Ai(r, t) = A4 core (r /lZ core (t)) a . From the 
definition of density this means: 

'^K^k)-' 3 f " r>a 

Now, if the density profile is to fall off as r increase as is the case in Nature, then {a — 3 < 0) 
which implies (a < 3). Combining this with (a > 0) we will have (0 < a < 3). Comparing this with 
the profile (p oc r~ a p), we have: (— a p — a — 3) and substituting this into (0 < a < 3), one obtains 
(0 < 3 — Up < 3). From (3 — a p < 3), we have (a p > 0), and from (0 < 3 — a p ), we have (a p < 3), 
hence (0 < a p < 3). 

Now, in this serene molecular, a small lamp begins to form - let this lamp have a radius lZi amp (t) 
and mass A4i am p- 1 shall pause a question: do we expect this lamp to cause any fundamental changes to 
the mass distribution Ai (r, t) = ar a + bl I think not. If this is the case, then our mass distribution must 
now be defined up till the radius of the lamp, i.e., A4(lZi am p, t) = M-iamp and this condition leads to: 

b = Mlamp - dT^Tamp^ henCe thuS: 

M(r,t) = a(r 3 - a r-n 3 l ^p+M larnp for r>K lamp (t). (26) 
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where we have substituted a = 3 — a p . Now inserting the condition that A4(7Z CO re,t) — M. core , we 
will have: 



M(r,t) - Mlamp 



'(*)-X7£ 



for r > Hi amp (t). (27) 



and putting all this together we will have: 



(3 — cxp TJ? Q p \ 
3 _ a ~ la ™* a \+ Miamp for r>K lamp (t). (28) 

TZcore" (t) - Ti-iamp J 

where Ai cs i(t) = M. core ~ ■Miampif)- Comparison of the above with (l24l i shows that the lamp in the 
above formula is the star in d24l i. 

We are certain the reader will have no problem with (|28T i because the lamp does not disrupt the mass 
distribution since it has no radiation hence we would expect a continuous distribution of mass right up 
the surface of the lamp as material will be flowing into the lamp. But this same lamp is a protostar and 
at somepoint it must switch on to become a star. At this moment, assuming the correctness of the thesis 
that at 8 — IOMq, the radiation field begins to push material away from the nascent star, we could from 
logic expect that the mass distribution must be continuous up till that time disruption starts at 8 — IQMq. 
During the time when the lamp's (or protostar's) mass is in the range < M.iamp{t) < 8 — 10.Mq, 
the MDF ( f28l ) must hold. From this, we have just justified the formula d24l i for the mass range: < 
■M s tar{t) < 8 — 10.M©. When the radiation field begins to be significant, we shall have to check and 
revise this formula. 

Now, the MDF d24l i. the gravitational field intensity - at any given time t and at any given point r 
inside the core from the surface of the star; will be given by: 

Circumstellar Gravitation „. , ~ .. .. 

Star s Gravitation 



= _ (GM l{ t)\ ( r^-nl^ { t) \ _ ,GM T (t)\ 

{ ' V n )\nlor7{t)-nlZ p (t)) V ^ ; 

Clearly, we have been able to separate the gravitation due to the star and the circumstellar material. 
Now, from the above, the inequality © becomes: 

GM csl (t) \ ( r^-^(t) \ ( GM star (t) \ Keff C atar (t) 

r 2 )\nlor a /{t)~nl; a a /(t)) \ n J a^c ' 

where the first term on the left hand-side of (f30T > is clearly the gravitational field intensity of the circum- 
stellar material and the second term is the gravitational field of the nascent star. 

5 RADIATION CAVITY 

The inequality (|5]l gives us the condition that must be met before the radiation field is powerful enough 
before it can push away (all) the circumstellar material inside the shell of radius r. Beyond this radius, 
the radiation field is not at all powerful enough to overcome the gravitational field. Unfortunately, one 
can not deduce this radius r from §5}. Fortunately, the inequality < f30b ; as does ©, (f30b tells us the 
conditions to be met before the radiation field is powerful enough to halt in-fall - in addition to this, 
( f30b sheds more information than (0 because in d30b we have quantified the MDF for the circumstellar 
material and this allows us to compute the region r where the radiation field is much stronger than the 
gravitational field. From ( f30b we deduce that the radiation field will create a cavity in the star forming 
core; in this cavity, the radiation field is much stronger than the gravitational field thus there will be 
found there-in no material but radiation hence the term - radiation cavity. To see this, that (f30b entails a 
cavity, we simple have to write (f30T > with r as the subject of the formula; so doing one arrives at: 



The Radiation Problem of Massive Stars Revisited 



11 




Accretion disk 



Fig. (1) For a non-spinning core at ~ 8 — 10M.Q, the nascent stars's accretion is halted (and 
importantly; in-fall is not reversed but only halted) because when the radiation field tries to 
create a cavity in which process the star is separated from its accretion source which is the 
circumstellar material; this means the star's mass accretion is halted because its mass can no 
longer grow since there exists no other channel(s) via which its mass feeds. Should the star's 
mass fall below ~ 8 — 10.M©, the circumstellar material will fall onto the nascent star until 
its mass is restored to its previous value of ~ 8 — 1QM.Q. In order for the radiation field to 
start pushing the circumstellar material, its mass must exceed ~ 8 — ICMq. Since there is 
no way to do this, in-fall is only halted and not reversed. Hence thus, the star's mass for a 
non-spinning star stays put at ~ 8 — lO-A/f©. As urged in Nyambuya j2010fet . this scenario is 
different for a spinning star because the ASGF (which comes about due the spin of the nascent 
star) allows matter to continue accreting via the equatorial disk inside the cavity as illustrated 
above. The accretion disk will exist inside the radiation cavity and this disk should according 
to the azimuthally symmetric theory of gravitation (Nyambuya 2010b), be channel mass right 
onto the nascent star right-up to the surface of the star without radiation hindrances. 



r > 



(KeffCstarit) - 4ircGM s tar(t)) (ill,"/ (t) 



TV. 



'(*) 



AncGMcsiit) 



n 3 s ;^(t) 



'^-cav [") 



(31) 

where TZ cav (t) is the radius of the cavity. Now that there is a cavity lets pause so that we can revise the 
MDF. Clearly, in the case where there are outflows, this must be given by: 



„3— a„ 



u 



3— a 



M(r, t) = M cs i(t) I -5— '" cu 3 "_ a + M.{t) for r > K cav {t). (32) 

\ f^-core V / '^-cav / 

where M cs i{t) = M core -M*(t) and M*(t) = M star (t) + Mdisk(t) + M utf(t): M dls k(t) is the 
disk mass inside the cavity at time t and M. utf the bipolar outflow contained in the cavity at time t. 
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Now, what this inequality ( f3TT > is "saying" is that, at any given moment in time when-after the star 
has surpassed the critical mass (8 — IOMq), there will exist a region r < H C av(t) where the radiation 
field will reverse the radially in-falling material and in the region r > TZ cav {t), for material therein, 
the radiation field has not reached a state where it exceeds the gravitational field hence in-fall reversal 
in that region has not been achieved. This region [i.e. r < TZ cav (t)] grows with time thus the radiation 
field slowly and gradually pushes away the material further and further away from the nascent star until 
TZ.cav{t) = TZ-ci where radial in-fall is completely halted and this will occur when the star has reached the 
critical core luminosity C* ore . The condition when the critical core luminosity has been shown earlier 
to lead to ( f]~2T > which is a Larson-like relation i.e., ( fT3l l, ipso facto, this strongly suggests that Larson's 
Law may not be a result of statistical sampling but a statement about and as-well a fossil record of the 
battle of forces between gravitation and the radiation field. 

By saying that the nascent massive star will create a cavity, we have made a tacit and fundamental 
assumption that its mass will continue to grow soon after the cavity begins to form and that its mass will 
thereafter continue to grow while in the cavity. But how can this be since the cavity separates the nascent 
star from the circumstellar matter? The nascent star is now without a channel to feed its mass so there 
can be no growth in its mass unless there exists a channel via which its mass feeds. At this juncture, we 
direct the reader to the readings Nyambuya (2010a, 20106). 

In Nyambuya <2010ab , as already said in the introductory section, we set-up the ASTG where-in 
the thesis was advanced to the effect (1) that, for a non-spinning star, its gravitational field is spherically 
symmetric (to be specific, it is only dependent on the radial distance from the central body); (2) that, for 
a spinning gravitating body, the gravitational field of that body in question is azimuthally symmetric, i.e., 
it is dependent on the radial distance (r) from the central body and as-well the azimuthal angle (6). In a 
follow-up reading i.e., Nyambuya (20106); we showed that the ASTG predicts (1) that bipolar outflows 
may very well be a purely gravitational phenomenon (i.e., a repulsive gravitational phenomenon) and 
also that; (2) along the spin-equator (define in there-in Nyambuya 1201061 1 of a spinning gravitating 
body, gravity will channel matter onto the spinning nascent star via the accretion disk (lieing along the 
spin-equator) thus allowing stars beyond the critical mass 8 — IO.M0, to come into being. It should be 
said that, accretion discs can also be formed by a number of different mechanisms other than the an 
Azimuthally Symmetric Gravitational Field (ASGF). 

The accretion of matter beyond the 8 — IOMq limit must only be possible for a spinning star 
because it possesses the ASGF that is needed to continue the channeling of matter onto the star via 
the accretion disk - see illustration in figure |(T)| For a non-spinning core the nascent stars's accretion 
can not proceed beyond 8 — IO.M0, it is halted because the moment the radiation field tries to create 
a cavity at the moment when the (non-spinning) star's mass is 8 — IOjW©, the (non-spinning) star 
in that very moment becomes separated from the surrounding circumstellar material. This means the 
(non-spinning) star's mass accretion is halted because its mass can no longer grow since there exists 
no other channel(s) via which its mass feeds. Should the (non-spinning) star's mass fall below 8 — 
IO.M0, the circumstellar material will fall onto the nascent (non-spinning) star until its mass is restored 
to its previous value of 8 — IOMq. This means the star's mass for a non-spinning star stays put at 
8 — IO.M0. As explained in the above paragraphs, this scenario is different for a spinning star because 
the ASGF (which comes about due the spin of the nascent star) allows matter to continue accreting 
via the equatorial disk. The accretion disk will exist inside the radiation cavity and this disk should 
according to the ASTG (Nyambuya 20106) channel mass right up to the surface of the star without 
radiation hindrances. The scenario just present is completely different from that projected in much of 
the wider literature where at 8 — 10 M.q, suddenly the radiation is so powerful it reverses any further 
in-fall. It is bona-fide knowledge that star formation is not a spherically symmetric process and from the 
above, it follows that stars beyond the 8 — lftM© limit must from with no hindrance form the radiation 
field and the only limit to their existence is gravitationally bound core with enough mass to form them. 
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6 DISCUSSION & CONCLUSIONS 

This contribution coupled with Nyambuya (20106) seem to strongly point to the possibility that the 
radiation problem of massive stars may not exist as previously thought. In the present reading, we find 
that beginning at the time when M. s tar{t) — 8 — lOA^©, the radiation field will create a cavity inside 
the star forming core and the circumstellar material inside the region TZ cav (t) < r < 1Z CO re (t) is going 
to be pushed gradually (importantly not blown away) as the radiation field from the star grows until a 
point is reached when the cavity is the size of the core itself, at which point complete in-fall reversal is 
attained. If the radiation field of the star is to grow, its mass must grow, thus, the cavity must not prevent 
accretion of mass onto the nascent star and this is possible for a spinning massive star. Once the cavity 
is created, the mass of the nascent will - for a spinning massive star; feed via the accretion disk and 
this disk is not affected by the radiation field. By saying the disk is not affected by the radiation field 
we mean the material on the disk is not going to be pushed away by the radiation field as it pushes the 
other material away because the azimuthally symmetric gravitational field of the star is powerful enough 
along this plane to overcome the radiation field - this has been shown or argued in Nyambuya (1201061) 
that this must be the case. 

The ASGF is only possible for a spinning star; since all known stars are spinning, every star should 
according to the ASTG have the potential to grow to higher masses. This means, massive stars should 
come into being because of their spin which bring about the much needed ASGF. A none spinning will 
will have no ASGF, hence no disk around it hence no channel via which to feed once the radiation field 
begins its toll. In this case of a none spinning star, this means once the star has reached the critical mass 
~ 8 — 1QM.Q, its mass can not grow any further because the moment it tries to grow, the star and the 
circumstellar material become separated due to the radiation field which in this case is stronger that the 
gravitational field. In this event, any further growth in mass of the star is annulled. This in actual fact 
means that for as long as there is circumstellar material, the mass of a none spinning star will stay put 
at ~ 8 — lOA^o because the moment it falls of slightly below ~ 8 — 10.M©, gravity becomes more 
powerful thus accrete mass only to restore it to its previous value of ~ 8 — lOA'f ©. In this case, we have 
an "eternal" stalemate between the gravitational and radiation field. 

An important and subtle difference between the present work and that of other researchers (Larson 
& StarfieldQUQ KahnHMl Yorke & KriigelUMl Wolfire & Cassinelli[TH571 Palla & Stahler[MJS 
Yorke [20021 Yorke & Sonnhalter 120031 ) is that we have seized on the observational fact that molecular 
clouds and cores are found exhibiting well defined density profiles. From this we computed the MDF 
which enabled us to find exactly the physical boundaries where the gravitational field is expected to 
be much stronger than the radiation field once the star exceed the critical mass. Additionally and more 
importantly is that from Nyambuya (20106) we have been able to argue that even after the cavity has 
been created mass will be channeled on to the star via the accretion disk. Without the ideas presented in 
Nyambuya ( 1201061 ). we would have been stuck because we where going to find ourself without a means 
to justify how the mass accretion continues once the cavity has been created. 

Importantly, we have pointed out a real problem in Yorke ( 12002b , Yorke & Sonnhalter (2003) and 
Zinnecker & Yorke (20071 , namely that these researchers have neglected the treatment of the circum- 
stellar material in their theoretical arguments leading to their definition of the radiation problem because 
they used Newton's inverse square law which clearly applies to a non-rotating mass in empty space - 
i.e., the inequality (0]i applies only for a star in empty space. In empty space, it is correct to say that the 
radiation field for a star of mass 10 M.q and beyond, will exceed the gravitational field everywhere in 
space beyond the nascent star's surface, but the same is not true for a star submerged in a pool of gas as 
is the case with the stars that we observe. 

Another important outcome is that it appears that Larson's Laws may well be a signature and fossil 
record of the battle of forces between the radiation and gravitational field. At present, it is thought of 
as being a result of statistical sampling thus the present brings us to start rethinking this view. We are 
not persuaded to think this is a result of statistical sampling. This view finds support from Weidner et 
al. (P2009)'s most recent and exciting work. In this work, these researchers present a thorough literature 
study of the most-massive star in several young star clusters in order to assess whether or not star 
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clusters are populated from the stellar initial mass function (IMF) by random sampling over the mass 
range (0.01.M© < M s tar < 150.Mo) without being constrained by the cluster mass. Their data 
reveal a partition of the sample into lowest mass objects (A4 c i < lOOA'fo), moderate mass clusters 
(100.M© < Md < IOOO.M0) and rich clusters above (M c i > lOOO-M©) where M c i is the mass of 
the molecular cloud. Their statistical tests of this data set reveal that the hypothesis of random sampling 
is highly unlikely thus strongly suggesting that there exists some well defined physical cause. 

In closing, allow us to say that, we do not claim to have solved the radiation problem but merely 
believe that what we have presented herein - together with the readings Nyambuya (201 0al [201061 ); 
is work that may very well be a significant step forward in the endeavor to resolving this massive star 
formation riddle. 
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